subroutine comput_flux1
use var
implicit none

    integer m,j,i
    
        do m=1,4
            call flux1_m(m)
            call flux1_p(m)
            do j=1,gridsy
                do i=0,gridsx
                    flux1(i,j,m)=fp1(i,j,m)+fm1(i,j,m)
                enddo                
            enddo            
        enddo

end subroutine